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We discuss the critical properties of the three-dimensional NJL model at nonzero 
' temperature. We show that the Z2-symmetric model undergoes a second order 

phase transition with 2d Ising exponents and its critical region is suppressed by 
■ a factor X/-J Nf. We also provide numerical evidence that the C/(l)-symmetric 

C^") ' model undergoes a BKT transition in accordance with the dimensional reduction 

' scenario. 

^ ■ Phase transitions in QCD at nonzero temperature have been studied 

| intensively over the last decade. However, since the problem of chiral sym- 

' metry breaking and its restoration is intrinsically non-perturbative, the 

number of available techniques is limited and most of our knowledge about 
this phenomenon comes from lattice simulations. 

Although there is little disagreement that the chiral phase transition 
in QCD with two massless quarks is second order, no quantitative work 
or simulations have been done that decisively determine its universality 
class. Universality arguments are very appealing due to their beauty and 
simplicity and they can be phrased as follows: At finite T phase transitions, 
the correlation length diverges in the transition region and the long range 
behavior is that of the (d— 1) classical spin model with the same symmetry, 
because the IR region of the system is dominated by the zero mode of 
the bosonic field. The contribution of non-zero modes does not affect the 
critical singularities but can be absorbed into non-universal aspects of the 
transition. Consequently, fermions which satisfy anti-periodic boundary 
conditions and do not have zero modes are expected to decouple from the 
scalar sector. A possible loophole to this standard scenario is that the 
mesons are composite // states and their size and density increase as T — > 
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T c . Therefore, if the transition region can be described as a system of highly 
overlapping composites the violation of the bosonic character of the mesons 
may be maximal and the fermions become essential degrees of freedom 
irrespective of how heavy they are. 

The NJL model has been proved to be an interesting and tractable 
laboratory to study phase transitions both numerically by means of lattice 
simulations and analytically in the form of large- Nf expansions 1,2 . The 
lagrangian density of the ?7(l)-symmetric model is 

C^^ + a + l^TT^ + ^i^+TT 2 ), (1) 

where the index i runs over Nf fermion flavors. There are several fea- 
tures which make this model interesting for the modelling of strong in- 
teractions: (i) The spectrum of excitations contains both "baryons" and 
mesons, namely the elementary fermions / and the composite // states; 
(ii) for sufficiently strong coupling g 2 > g 2 it exhibits spontaneous chiral 
symmetry breaking implying dynamical generation of a fermion mass; (iii) 
for 2 < d < 4 there is an interacting continuum limit at a critical coupling, 
which for d = 3 has a numerical value g 2 w 1.0 /a in the large- Nf limit if 
a lattice regularization is employed. At leading order in 1/Nf the model 
undergoes a second order transition at T c = 3 (Mf is the dynamical 
fermion mass at T = 0) with Landau- Ginzburg mean field scaling. This 
result is valid for both continuum and discrete chiral symmetries. 

We studied the critical behavior of the ^-symmetric model with finite 
Nf by performing lattice simulations near T c 4 . We simulated the model 
exactly for Nf — 12 with the Hybrid Monte Carlo method. The staggered 
fermion lattice action and further details concerning the algorithm can be 
found in 2 . The temporal lattice size was L t — 6 and the spatial size varied 
from L s = 18 to 50. The expectation value of the auxiliary sigma field serves 
as a convenient order parameter for the theory's critical point. By using 
finite size scaling (FSS) techniques we extracted the critical exponents. The 
results which are summarized in Table 1 support the dimensional reduction 
scenario, because the values of the exponents are in good agreement with 
those of the 2d Ising model rather than the mean field theory ones. This 
implies that the role of composite mesons is relevant near T c even when Nf 
is large but finite. Simulations with the meron-cluster algorithm of the 4c? 
NJL model gave 3d Ising exponents with very high accuracy. 5 . 

Next we tried to understand how the large- Nf mean field theory predic- 
tion reconciles with the dimensional reduction and universality arguments. 
The answer is that the large- Nf description has its applicability region. We 
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studied this in detail for a Yukawa theory 4 and we showed that the phe- 
nomenon which leads to an apparent contradiction is the suppression of the 
width of the non-mean field critical region by a power of 1/Nf. The 2d Ising 
critical behavior sets in when T ^> M a (T) (M a (T) is the thermal mass of 
the a meson). If the cutoff A ^> T the renormalized self-interaction cou- 
pling A(T) in the large- Nf limit is close to the IR fixed point of the Yukawa 
theory and is given by A(T) <~ T 4 ~ d /Nf for 2 < d < 4. The mean field 
approximation breaks down because of self-inconsistency when the value of 
the coupling of the (d—1) dimensional scalar theory \d-i ~ TX(T) on the 
scale M^~ d (T) (the power d — 5 comes from dimensional analysis) is not 
small anymore. Therefore, for d = 3 the Ginzburg criterion for the applica- 
bility of the mean field scaling is given by M a {T) 3> Tj JNf. This scenario 
was also verified numerically by measuring (a) = Mf(T) = M a (T)/2 at the 
crossover into the mean field region 4 . 

Additional evidence in favor of the dimensional reduction scenario was 
produced in studies of the J7(l)-symmetric NJL3 model 6 ' 7 - 8 . Both analyt- 
ical and numerical results showed that its phase structure at T ^ is the 
same as that of the 2c? XY model. As we already mentioned, at leading 
order in 1/Nf this model undergoes a second order chiral phase transi- 
tion . This conclusion, however, is expected to be valid only when Nf is 
strictly infinite, i.e. when the fluctuations of the bosonic fields are neglected, 
otherwise it runs foul of the Coleman-Mermin- Wagner (CMW) theorem, 
which states that in two-dimensional systems the continuous chiral symme- 
try must be manifest for all T > 0. Ncxt-to-leading order calculations 6 ' 8 
demonstrated that the model undergoes a Berezinskii-Kosterlitz-Thouless 
(BKT) transition 10 , in accordance with the dimensional reduction scenario. 
It is easy to visualize the BKT scenario if we use the "modulus-phase" 
parametrization: a + in = pe ld . In two spatial dimensions logarithmically 
divergent infrared fluctuations do not allow the phase 9 to take a fixed di- 
rection and therefore prevent spontaneous symmetry breaking via (9) 7^ 0. 
The critical temperature Tbkt is expected to separate two different chi- 
rally symmetric phases: a low T phase, which is characterized by power 
law phase correlations ( e te ^e~ te ^) ~ x~' n( - T ^ at distances x > 1/T but 
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no long range order (i.e. a spinwave phase where chiral symmetry is "al- 
most but not quite broken"), and a high T phase which is characterized 
by exponentially decaying phase correlations with no long range order. In 
other words for T < Tbkt there is a line of critical points characterised by 
r](T) = ^ 9 forT« Tbkt and at T = T BKT v(T) is expected to be 0.25. 

In order to study the behavior of the chiral symmetry at T > in the 
absence of a fermion bare mass in the Lagrangian the best thing to measure 
is an effective "order parameter" |$| = V 'a 2 + 7r 2 , which is a projection onto 
the direction of <& Q = (a, w) separately for each configuration. In Fig.l we 
plot |$| versus (5=1/ g 2 for fixed lattice temporal size L t — 4 and different 
lattice spatial sizes together with the results from simulations of the model 
with a Z2 chiral symmetry. In both cases Nf = 4. At the bulk critical 
coupling /^ ulk = l/g 2 «0.86 7 the lattice spacing becomes zero and T^oo. 
It is clear that the order parameter of the Z 2 model is independent of the 
lattice size until just before the transition at /3,f 2 =0.565(3), whereas in the 
U(l) model |$| has a strong size dependence for a large range of values of 
/?, i.e. it decreases rapidly as the spatial volume increases in accordance 
with the expectation that chiral symmetry should be restored for T > 0. 
The finite spatial extent L s provides a cut-off for the divergent correlation 
length and according to the BKT scenario the slow decay of the correlation 
function (e l9 ^ e~ ie ^) for T < Tbkt ensures a non-zero magnetization 
even in a system with very large size. In Fig.l we plot the susceptibility 
of the order parameter \ = ^((l^l 2 ) — (l^l) 2 ); measured on lattices with 
different spatial size. We observe from this figure that: (a) there is a phase 
transition since the peak of x diverges as we increase L s ; and (b) in the low 
T phase the susceptibility has a stronger size dependence than in the high 
T phase and the errors at low T are much larger than at high T. This is 
strong evidence that the system is critical in the low T phase in accordance 
with the BKT scenario. It is instructive to compare the results for the U(l) 
model with results extracted from recent simulations of the SU{2) x SU(2)- 
symmetric N JL3 model 12 . In the latter case the correlation length is finite 
for T > and the effective order parameter goes to zero more rapidly as 
we increase L s than in the U(l) case. It was also shown that both models 
have a pseudogap phase where the fermions remain massive in the absence 
of long range correlations 6 > 7 ' 12 . Recently, simulations of an extension of 
the attractive Hubbard model with the meron-cluster algorithm confirmed 
with high accuracy that < i] < 0.25 11 . 

In summary, we have shown that the dimensional reduction scenario 
describes correctly the phase transition of simple fermionic field theories 
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Figure 1. Order parameter and susceptibility vs. coupling fi for various values of L. 



with composite mesons. We plan to extend our work to the four-dimensional 
NJL model and study the dependence of the width of the critical region 
on the degree of chiral symmetry. It has already been shown that in the 
(3 + l)d O(N) sigma model the critical region is much wider for N > 2 
than for the discrete symmetry case 13 . Finally, it will be interesting to 
understand the issue of universality in two-flavor QCD. Simulations with 
massless quarks of QCD with a four-fermi interaction 1 near the continuum 
limit may provide decisive results in this direction. 
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